Negative Refraction Requires Strong Inhomogeneity 
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The paper establishes explicit lower bounds for the lattice cell size of periodic structures (meta- 
materials and photonic crystals) capable of supporting backward waves and producing negative 
refraction. At optical frequencies, this result implies strong inhomogeneity, in the sense that the 
cell size cannot be negligible relative to the vacuum wavelength and the Bloch wavelength. 



PACS numbers: 41.20.-q, 41.20.Cv, 41.20.Jb, 42.25.Bs, 42.70.Qs 



■ 

o : 
o . 

(N ■ 
CD ■ 

m ; 
oo ■ 

<n : 



. 

o . 
^— > ' 

Or 

q : 
o ■ 

CZJ ■ 

>v 
^ : 



> 



o 
o 



o 
o 



X 



Negative refraction (electromagnetic waves bending 
the 'wrong' way at material interfaces) and the closely 
related phenomenon of backward waves (phase velocity 
at an obtuse angle with group velocity) have become one 
of the most intriguing areas of research in nanophotonics 
this century, with a number of books and review papers 
readily available, e.g. [J H S3, S3, S3, S3, and hun- 
dreds of research papers published. As early as in the 
1940s, Mandelshtam pointed out [13] that waves would 
refract negatively at the interface boundary between a 
regular and a backward-wave medium. In 1967, Veselago 
showed that media with simultaneously negative (rela- 
tive) dielectric permittivity e r and magnetic permeabil- 
ity [i r would support backward waves and exhibit other 
unusual behavior of wave propagation and refraction [33] • 
In 1999-2000, Pendry et al. [Il| proved theoretically 



and Smith et al. demonstrated experimentally nega- 
tive refraction in an artificial medium with split-ring res- 
onators. Furthermore, Pendry discovered that Veselago's 
unusual 'lens' - a slab of a negatively refracting material 
- could produce a perfect image of a point source, thereby 
beating the diffraction limit [281 ] . 

Truly homogeneous materials, in the Veselago sense, 
are not currently known. Consequently, much effort has 
been devoted to the development of artificial metama- 
terials capable of supporting backward waves and pro- 
ducing negative refraction S !, 0, S S3, 0, S3 SI Si . 

Separately from the progress in metamaterials, negative 
refraction has been observed and analyzed in singly and 
doubly periodic waveguides |38| an d in photonic crystals 
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All these intriguing findings have led to the presump- 
tion that there are two species of negative refraction, one 
occurring in photonic crystals and another one in meta- 
materials. Conceptually, the latter are viewed as proto- 
typical 'Veselago media'. 

There are, indeed, salient differences between metama- 
terials and crystals in terms of the underlying structure, 
composition and fabrication (e.g. lossless dielectric inclu- 
sions vs. lossy metallic resonators of various kinds). On 
a more fundamental level, however, all such structures 
can be characterized by a periodically varying complex 
dielectric function, and from that point of view it is legit- 
imate to examine possible principal differences between 



metamaterials and photonic crystals. Importantly, can 
metamaterials, as a matter of principle, be (arbitrarily) 
close to an ideal homogeneous Veselago medium? 

In the experimental and computational examples of 
the references cited above, the cell size as a fraction of 
the vacuum wavelength varies between ~ 0.11 -j- 0.42. 
One would hope that further improvements in nanofab- 
rication and design could bring the cell size down to a 
smaller fraction of the wavelength, thereby approaching 
the Veselago case of a homogeneous material. However, 
the main conclusion of the present paper is that the cell 
size is constrained not only by the fabrication technolo- 
gies but by fundamental lower bounds as well. 

The analysis in this paper relies on the usual 2D and 
3D renditions of time- harmonic Maxwell's equations, and 
it is assumed that bulk material parameters are applica- 
ble with a reasonable level of accuracy. At optical fre- 
quencies, the relative intrinsic permeability of all media 
can be set to unity ([13], §60; [13].) To streamline the 
mathematical development, we focus on square / cubic 
Bravais lattice cells with size a in 2D/3D and introduce 
dimensionless coordinates x = x/a, etc., so that in these 
tilde-coordinates the 2D / 3D problem is set up in the 
unit square / cube. The s-mode in the tilde-coordinates 
is described by the 2D wave equation 
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where E is a one-component electric field phasor and 
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Here c and Ao are the speed of light and the wavelength 
in free space, respectively. The relative permittivity e r 
is a periodic function of coordinates over the lattice. In 
3D, the governing equation is 



V x V x E 



t!> 2 eE 



(3) 



The fundamental solutions of the field equation in pe- 
riodic structures are known to be Bloch-Floquet waves 
with a (yet undetermined) Bloch vector K b '■ 



E(r) = EpER(f) exp(iK B • f ) 



(4) 



where f is the position vector. Subscript 'PER' im- 
plies periodicity with respect to any lattice vector 
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(n x a,n v a,n z a) with integer n x , n y , n z in the 3D case. 
Eper can be expanded into a Fourier series 

EpER(f) = y~] e n exp(i27rn • f), (5) 

where e n are the Fourier coefficients and index n runs 
over the integer lattice Z 2 or Z 3 in 2D/3D. 

For analysis and physical interpretation of energy flow, 
phase velocity and other properties of the Bloch wave, 
it is convenient to view it as a suite of spatial Fourier 
harmonics (plane waves) [3J. From (01 and ([5]), 



E(f) 



E n 



e n exp(i27rn- f ) exp(iKs • f) (6) 



The decomposition of the magnetic field is similar. 

It is important to note from the outset 0] that the indi- 
vidual plane-wave components E n of the electromagnetic 
Bloch wave do not satisfy Maxwell's equations in the pe- 
riodic medium and therefore do not represent physical 
fields. Only taken together do these Fourier harmonics 
form a valid electromagnetic field. 

It is straightforward to verify that the plane waves in 
the decomposition are orthogonal functions over the lat- 
tice cell (in the sense of standard vector L2 inner prod- 
uct). Hence, by Parseval's theorem, the time- and cell- 
averaged Poynting vector <P> = i<Rc{ExH*}> 
can be represented as the sum of the Poynting vectors 
for the individual plane waves Q: 



< P > = 



\e n \ 2 



(7) 



Group velocity doj/dk is clearly the same for all plane 
wave components, and hence group velocity for the whole 
Bloch wave can be defined as v g — Ooj/OKb- In cases 
of weak dispersion, this velocity indeed approxim ately 
represents signal velocity in the periodic medium [34ll37|. 

It is well known that in Fourier space the scalar wave 
equation (J]) becomes 



IK, 



27rn| 2 e„ = Cj 2 e n _ 



(8) 



where e n are the Fourier coefficients of the dielectric per- 
mittivity e: 



Ee n exp (i27rn • f ) 
1 1 



(9) 



Indeed, the right hand side of © is Fourier-space con- 
volution corresponding to real-space multiplication eE. 
The left hand side represents —V 2 . 

Refraction at the interface between the periodic struc- 
ture and air (or another homogeneous dielectric) depends 
not only on the intrinsic characteristics of the Bloch wave 
in the bulk, but also on the extrinsic conditions at the 
interface boundary - namely, the 'excitation channel' 
i.e. the Fourier component of the Bloch wave that cou- 
ples to the incident wave in the air 0, 13, 15|. Intrinsic 
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FIG. 1: The photonic band diagram of the Gajic el al crystal. 
TE modes (p-polarization, one-component H field) - squares, 
solid lines. TM modes (s-polarization, one-component E 
field) - circles, dashed lines. 



properties include the forward or backward character of 
the wave - that is, whether the Poynting vector and phase 
velocity (if the latter can be properly defined) are at an 
acute or obtuse angle. 

For illustration and further analysis, it is convenient to 
have a specific example in mind (however, the analysis 
and conclusions will be general). Consider the structure 
proposed by R. Gajic, R. Meisels et al [TH, [l6j]. Their 
photonic crystal is a 2D square lattice of alumina rods 
in air. The radius of the rod is r ro( j = 0.61 mm, the lat- 
tice constant a = 1.86 mm, so that r ro d/a ks 0.33. The 
band diagram, computed using the plane wave method 
with 441 waves for s- and p-modes appears in Fig. [TJand, 
apart from the scaling factors, is very close to the one in 
TBI 16| , where various cases of wave propagation and re- 
fraction are studied. In the context of this paper, of most 
interest is negative refraction for small Bloch numbers in 
the second band of the p-mode. 

We observe that the TE2 dispersion curve is mildly 
convex around the T point (Kb = 0, uja/2irc w 0.427), 
indicating a negative group velocity for small positive Kb 
and a possible backward wave. 

An additional condition for a backward wave must also 
be satisfied: the plane-wave component corresponding to 
the small positive Bloch number must be appreciable (or 
better yet, dominant). The distribution of Poynting com- 
ponents of the same wave is shown in Fig. [21 It is clear 
from the figure that the negative components outweigh 
the positive ones, so power flows in the negative direction. 
(Details can be found in [13|.) 

However, the normalized band diagram indicates that 
negative refraction disappears in the homogenization limit 
when the size of the lattice cells tends to zero, provided 
that other physical parameters, including frequency, are 
fixed. Indeed, the homogenization limit is obtained by 
considering the small cell size - long wavelength limit 
a — > 0, K — > (see 0, H2| for additional mathemati- 
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FIG. 2: The plane- wave Poynting components P m for the 
Gajic et al crystal (arb. units). Second ff-mode (TE2) near 
the r point on the V — > X line. 



cal details on Floquet-based homogenization theory for 
Maxwell's equations). As these limits are taken, the 
problem and the dispersion curves in the normalized co- 
ordinates remain unchanged, but the operating point 
(cj,K) approaches the origin along a fixed dispersion 
curve - the acoustic branch. In this case phase veloc- 
ity in any given direction I, lo/Ki = Co/ Ki, is well defined 
and equal to group velocity duo / dKi simply by definition 
of the derivative. No backward waves can be supported 
in this regime. 

This conclusion is not surprising from the physical per- 
spective. As the size of the lattice cell diminishes, the 
operating frequency increases, so that it is not the ab- 
solute frequency uo but the normalized quantity Co that 
remains (approximately) constant. Indeed, a principal 
component of metamaterials with negative refraction is 
a resonating element 0, M, ED, [33l whose resonance fre- 
quency is approximately inverse proportional to size . 

It is pivotal in this paper to make a distinction be- 
tween strongly and weakly inhomogeneous cases of wave 
propagation. The latter is intended to resemble an ideal 
'Veselago medium,' with the Bloch wave being as close 
as possible to a long-length plane wave. Toward this end, 
the following conditions characterizing the weakly inho- 
mogeneous backward-wave regime are put forth: 



• The first-Brillouin-zone component of the Bloch 
wave must be dominant; this component then de- 
fines the phase velocity of the Bloch wave. 

• The other plane-wave components collectively pro- 
duce energy flow at an obtuse angle to phase veloc- 
ity. 



• The lattice cell size a is small relative to the vacuum 
wavelength Ao ; o/Xq <C 1 . 

• At the air-material interface, it is the long- 
wavelength, first-Brillouin-zone, plane wave com- 
ponent that serves as the excitation channel for the 
Bloch wave. 

If any of the above conditions are violated, the regime 
will be characterized as strongly inhomogeneous: the 
EM wave can "see" the inhomogeneities of the mate- 
rial. By this definition, in the weakly inhomogeneous case 
the normalized Bloch wavenumber K B must be small, 
K B = K B a <C ir. Larger values of K B would indicate 
a strongly inhomogeneous (or, synonymously, 'photonic 
crystal' or 'grating') regime, where the lattice size is com- 
parable with the Bloch wavelength. As we shall see, 
under reasonable physical assumptions, backward waves 
cannot be supported in the weakly inhomogeneous case; 
strong inhomogeneity is required. 

As a preliminary step in the analysis, it is instructive 
to examine the direction of power flow for small Kb in 
the lossless cae (real e). The average Poynting vector is, 
according to ((7]) and with a convenient normalization, 



P = 2u)/i < P > = K E 



e {K L 



m—1 



{K B +2imi) 



- (K B -2nm) e~ m {K B ) 
(10) 

The scalar form is used for notational convenience only; 
the vectorial case is quite similar. It is, however, essential 
to indicate explicitly that the Fourier amplitudes e m de- 
pend on the Bloch parameter K B . Since the waves corre- 
sponding to ±K B are complex conjugates of one another, 
we have e- m (K B ) = e^—Ks), and the expression for 
the Poynting vector becomes 



P = K F 



e (K B ) +> e m (K B ) + e m {-K B 

£ — 'm— 1 \ 



e m {K B ) - e m (-K B ) 



(11) 



The first two terms in (llip are directly proportional to 
K B . To make this small parameter explicit in the third 
sum as well, we write 



P = K f 



leol 2 



Eoo 
m=l 



27T 



e m (K B ] 



dK E 



e m (-K B ) 



(12) 



For small Co, the positive term j eo | 2 in the square brackets 
tends to be dominant, making it difficult to produce a 
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negative power flow and a backward wave. This is so 
because the magnitudes of all spatial harmonics except 
for io are for small u) constrained by © : 



< 



CL> 2 



K, 



2nn 



= 1. 



(13) 

The arguments above suggest that there must be a 
lower bound for the relative cell size a/Xo = (D/27T when 
the medium could still support backward waves. To the 
best of my knowledge, this question has not so far been 
posed explicitly in the literature. 

In the remainder, we investigate the constraints on 
the periodic in the weakly inhomogeneous backward-wave 
regime. This implies that Kb — Kbcl C 1. To sim- 
plify mathematical analysis, we focus on the limiting case 
Kb = 0, but the conclusions will apply by physical conti- 
nuity, to small Kb ■ We first turn to the s-mode governed 
by the 2D equation (P). For Cj ^ and r\ 



~,-2 



u> 



eE = -fyV 2 £ 



(14) 



Further analysis relies on the inversion of V 2 . To do 
this unambiguously, let us split E up into the zero-mean 
term E±_ and the remaining constant E$: E — Eq + 
E±. Symbol '_L' indicates orthogonality to the null space 
of the Laplacian (i.e. to constants). To eliminate the 
constant component Eq , we integrate (fT4]) over the lattice 
cell. Integrating by parts and noting that the boundary 
term vanishes due to the periodic boundary conditions 
(K B = 0), we get 



En — —in 



eE ± ain, en^Q 



(The exceptional case eo = is mathematically quite 
intricate and may consitute a special topic for future re- 
search.) With Eq eliminated, the eigenvalue problem for 
E± becomes 



eE, dQ, 



= -r 1 V 2 E 1 



Since E±_ by definition is zero-mean, 

VI 2 {e [E ± - i^(e, E ± )]} = - V E A 



(15) 



where Vj_ is the zero-mean inverse of the Laplacian. 
Fourier analysis easily shows that this inverse is bounded 



< 



(the Poincare inequality): 
taking the norm of both sides of (fT"5|) , we get 

max V x I | c |max / leol) 



(4^ 2 



Then, 



(16) 



This result, that can be viewed as a generalization of the 
Poincare inequality to cases with variable e r , leads to a 
simple lower bound for the lattice cell size, with the mean 
and maximum values of e as parameters: 



a 



4tt 2 



> 



1 



</k~ |) 



(17) 



Turning now to the vector field formulation ([3]), we 
deal with 2D and 3D cases simultaneously and rewrite 
the field equation as 



eE = -r? V x V x E 



(18) 



The V x V x operator can be inverted unambiguously if 
the result, denoted with (Vx)J 2 , is sought in the func- 
tional space Hj_(f2) of divergence- free zero- mean fields. 
For any such field u, V x V x u = — V 2 u and hence 



V 2 u 



= (V x V x u, V x V x u) 



= (V 2 u, V 2 u) > (4tt 2 ) 2 (u, u) = (4tt 2 ) 2 ||u|| 2 

This implies that the inverse curl-curl, considered as an 
operator with its range in H^, is bounded: 



(Vx); 2 < (in 2 )- 1 



(19) 



The relevant splitting of E is into the zero-mean 
divergence- free term Ej^ 6 and the curl-free remain- 
der [40j Eo = — V<^o (the Helmholtz decomposition): 
E = ~E± — V<j)Q. Field E^ is in fact, up to the fac- 
tor iCj, the magnetic vector potential with the Coulomb 
(zero-divergence) gauge. 

Taking divergence (in the distributional sense) of the 
governing equation (|18p and integrating over the cell, one 
eliminates the electrostatic term V(f>o and arrives at an 
eigenvalue problem for Ej_: 



E 



v/r 1 v 



(eEj 



-7?V X V X Ej_ 



assuming that the electrostatic operator C t ~ V • eV 
is nonsingular. Equivalently, since E^ is by definition 
divergence-free and zero-mean, 



(Vx): 2 { e [E ± - V£r x V.(eE ± )]} 



VV± (20) 



an upper bound for r\ can be obtained by taking the L2- 
norms of both sides, with (1191) in mind: 



M < (4 7 r 2 )- 1 |e| max (l + |A| n , ax (£7 1 )|e| max ) (21) 

This estimate is analogous to the scalar one (TTB1 . ex- 
cept that the maximum eigenvalue |A| max (£7 1 ) ( n °t to 
be confused with the wavelength) appears instead of the 
inverse mean value |eo| _1 of the permittivity. This eigen- 
value is bounded unless the operating frequency is close 
to the quasi-static plasmon resonance value. In the most 
general situation, no simple estimate of |A| max ('C ( r 1 ) i s 
available, but it can be computed numerically using a 
number of algorithms (see e.g. for any given distri- 
bution of e in the lattice cell. 

At the same time, there are practically important situ- 
ations where the bound for rj can be made more explicit. 
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FIG. 3: Bounds on the normalized cell size and a few repre- 
sentative data points from the literature. 



One such case is that of non-plasmonic materials, when 
£ > £min > throughout the lattice cell. Then 

IAUxCO < e m in |A| maX (V: 2 ) < (47T 2 £min )- 1 

and from the estimate pTj) for the following bound on 
the normalized cell size emerges: 



4ir 2 \r]\ 



> 



ax (1 + |e|max/(47r 2 e mi „)) 
(22) 

The lower bounds (| 1 T|) and ([22]) are plotted as a func- 
tion of |e| m ax in Fig. [3j for e m j n = 1 and two values of 
eo (0.5 and 0.1). For illustration, several representative 
data points (both theoretical and experimental) from the 
literature are also shown in the figure. In the microwave 
regime, when metals are very good conductors and con- 
sequently |e| max is high, the theoretical bound for the 
cell size is non-restrictive and the respective data points 
(Smith, Shelby, Houck, and others) easily turn up above 
the relevant theoretical curve; these points lie off the 
chart in Fig. [3] 

The Cubukcu et al data point lies below the theoreti- 
cal line (j2"2"|) ; however, there is no contradiction because 
in this instance negative refraction occurs in the vicin- 
ity of the M point, where the Bloch wavelength and the 
lattice cell size are comparable. This constitutes, by our 
definition, a strongly inhomogeneous case to which the 
theoretical bound does not apply. 

The Moussa and Gajic data points for non-metallic 
crystals lie only slightly above the theoretical bound, 
indicating that this bound can be approached in some 
cases. Still, it must be stressed that the theoretical lim- 
its on the cell size are necessary, but in general not suf- 
ficient, conditions for negative refraction. A sufficiently 
large lattice cell size makes it possible for higher-order 
Fourier harmonics of the Bloch wave to outweigh the 



first-Brilloin-zone harmonic, but does not guarantee that 
they will do so and that they will have the desirable sign. 

Another case where the theoretical bound (f2"T|) can be 
made more explicit is that of a lossless host medium, 
e = eh, with embedded 'inclusions' e = ti = e- + ie" 
(spheres, split-ring resonators, fishnets, horseshoes, rods, 



etc.). The eigenvalue |A| max (£ e )~ 



I ALL (^e) can be 



estimated from the electrostatic energy functional 

(eV0, V0) - e h W h + (ej + ieJ')W< 

where W h ,i = J Qh . |V0| 2 dfl Then |(eV0, V(/>)| 2 is a 
quadratic form with respect to Wh,i and can be bounded 
by direct evaluation of its minimum eigenvalue. The end 
result, for small losses e" <C e;, + leA, is 



|A|min(£<0 



> 



2 //2 
t h t i 



2(^ + 



This estimate can be used in conjunction with the general 
bound ([21]). 

In summary, it has been proved that periodic struc- 
tures capable of supporting backward waves and pro- 
ducing negative refraction in the optical range must be 
strongly inhomogeneous. More precisely, the lattice cell 
size, as a fraction of the vacuum wavelength and/or the 
Bloch wavelength, must be above certain thresholds es- 
tablished in this paper. These thresholds contain the 
maximum, minimum and mean values of the complex di- 
electric permittivity as key parameters. In the presence 
of good conductors (e.g. at microwave frequencies) such 
theoretical constraints are not very restrictive. However, 
at optical frequencies and/or for non- metallic structures 
the bounds on the cell size must be honored and may 
help to design metamaterials and photonic crystals with 
desired optical properties. 
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